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Abstract We study the radial- velocity and light curves 
of the two eclipsing binaries EE Aqr and Z Vul. Us- 
ing the latest version of the Wilson & Van Hamme 
(2003) model, absolute parameters for the systems are 
determined. We find that EE Aqr and Z Vul are near- 
contact and semi-detached systems, respectively. The 
primary component of EE Aqr fills about 96% of its 
'Roche lobe', while its secondary one appears close to 
completely filling this limiting volume. In a similar 
way we find fill-out proportions of about 72 and 100% 
of these volumes for the primary and secondary com- 
ponents of Z Vul respectively. We compare our results 
with those of previous authors. 

Keywords Variable stars- Binaries- Eclipsing binary 
1 Introduction 

Eclipsing binary stars studies often involve the combi- 
nation of photometric (light curve) and spectroscopic 
(mainly radial velocity curve) data. The analysis of 
the light-velocity curves enable astronomers to obtain 
absolute physical parameters describing the system and 
its components. The physical parameters derived from 
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the photometric and spectroscopic data can help to im- 
prove our understanding of physical processes in stars. 
In view of the importance of EB-EA type binary stars 
in the problem of mass exchange in binaries and in the 
theory of their evolution, the study of these systems 
play an important role in finding of the complex initial 
phases of stellar evolution. In order to understand the 
physical nature of EB and EA type binaries, two can- 
didates, EE Aqr and Z Vul, were selected based on 
the following criteria: a) the shape of the light curve, 
b) their photoelectric measurements are not enough and 
also somewhat the analysis results derived by the previ- 
ous authors appear to differ from each other. The first, 
we state their history in literature and then describe 
our solution method. 

1.1 EE Aqr 

The binary system EE Aqr (HD213863, BD-20°6454, 
SAO 191236, 02000 = 22''34'"4P.87, ^2000 = 19°51'34") 
was discovered to be variable by Strohmeier & Knigge 
(1960). Strohmeier et al. (1962) derived its echpsing 
behavior with a orbital period of 0.51 d and classified 
it as an Algol (EA) type eclipsing binary. Later, the 
photoelectric observations were obtained by Williamon 
(1974), Padaha (1979) and Covino et al. (1988). WiUia- 
mon (1974) who showed the light curves to be more like 
that of 13 Lyrae (EB) type with a period of 0.5089951 d 
and analyzed these using the method of Russell & Mer- 
rill (1952). Thus preliminary photometric orbital pa- 
rameters of the system were derived. Slightly different 
orbital parameters were determined by Padalia (1979) 
for the system using the same method. Russo & Sol- 
lazzo (1982) found many inconsistencies in applying the 
Russell and Merrill (1952) method of solution. There- 
fore they reanalyzed both V light curves of Williamon 
(1974) and Padaha (1979) using the Wilson & Devin- 
ney (1971) model. They obtained the same results for 
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both light curves while the Russell-Merrill model does 
not. They found a solution as a semi-detached system. 
The system has also been studied spectroscopically by 
Hilditch & King (1988), from which radial-velocities 
were measured for both components by means of the 
cross-correlation code VCROSS. The first simultaneous 
solution of photometric and radial-velocity curves was 
performed by Covino et al. (1990) based on the Wil- 
son & Devinney (1971) model. Their solution indicated 
that EE Aqr might be a contact configuration not yet 
in thermal equilibrium. From 1990 to now, no analysis 
of light curves and new photometric observations has 
been reported. 

The period variations of the system has studied by 
Srivastava (1987) who found irregular changes of the 
order of ±10~^ d and many observations of the times 
of minima (Mallama (1980), Covino et al. (1988), Dccg 
et al.(2003)) so far. 

The spectral classes of the components of EE Aqr 
were given as FO (Williamon (1974)), FV-^A (Padaha 
(1979)), A8V-|-(K3-K4)(Russo & SoUazzo (1982)). 

1.2 Z Vul 

The variability of Z Vul (HD181987, aaooo = 19'' 21" 
39M1, (520 00 = 25°34'29".45) was discovered to be an 
eclipsing binary by Herschel (It has been reported by 
Astbury, 1909). Plaskctt (1920), Petrie (1950), Roman 
(1956), Popper (1957b) and Cester et al. (1977) de- 
termined the spectral type of the components as B3, 
B4+B6, B5V+A, B3-4V+A2-3III and B2V+A1III re- 
spectively. The first researches on the system presented 
its eclipsing behavior with a orbital period slightly less 
than 2.5 d and classified it as an Algol (EA) type eclips- 
ing binary. Cester et al. (1977) analyzed Broglia's 
photometric observations (1964) using the WINK pro- 
gram (Wood, 1972) and determined slightly different 
masses, radii and luminosities for the system with those 
found in the literature. They found a solution as a semi- 
detached system. 

Peters (1994) presented the first lUE observations to 
investigate the effect of the cross-section and tempera- 
ture of the primary on the mass transfer and mass loss 
in the system. They found that the primary involve 
with weak winds, originating from the gas stream will 
strike the primary instead of forming an accretion disk. 
The lack of H-alpha emission in system supports this 
proposition. 

The first results from far-UV observations obtained 

by Peters & Polidan (1997) about the nature of the cir- 
cum stellar material in the system. Their observations 
showed no infalling from a gas stream. So it is consis- 
tent with this fact that the system is sufficiently close 
and can not establish a disk. 



It is recognized that depend on applying present 
models, we can obtained the same or diS'erent results 
of analyzing light-velocity curves and thus it follows an 
low or high understanding of stellar evolution. Therefor 
it is need to research a tool for a better comprehension 
of the binaries structure and evolutive stage. 

In this study the authors adopted the more realistic 
close binary model based on the Roche model. In order 
to obtain accurate solutions in this situation, the light 
curves were analyzed by means of the latest the Wil- 
son's computer code. The computing code used was de- 
veloped by the Wilson & Van Hamme code (2003) (here 
after WV) to determine the parameters of the system. 
We selected it as our analysis research tool both for 
its intrinsic virtues and because of some improvements 
in comparison with earlier code (Wilson & Divinney, 
1971). The fourth revision (2003) is improved for ex- 
ample based on Kurucz's new atmospheres, log g as 
a parameter (allowing for handling giants, sub-giants, 
etc., in addition to main sequence stars) so that tem- 
perature ranges vary according to log g together with 
19 abimdances (relative to the sim). Hence we felt it 
useful to reanalyze these systems and get the most re- 
liable elements with the model of WV, which was de- 
veloped especially to handel close systems and compare 
our finding with results of the previous authors. 

All the data we analyze in the paper are taken 
from the literature so that these data appeared to be 
free from any peculiarities reported by the previous 
researches. We consider photometric observations of 
Williamon (1974) (here after WI) and spectroscopic 
observations of Hilditch & King (1988) (here after HK) 
for EE Aqr and the photometric observations of Broglia 
(1964) (here after BR) and spectroscopic observations 
of Popper (1957b) (here after PO) for Z Vul. 

The present paper is organized as follows. The as- 
sumptions are described in section 2. Photometric solu- 
tions of light curves is given in section 3. Spectroscopic 
solutions of radial velocities curves in section 4. Ab- 
solute elements for the primary and secondary compo- 
nents in section 5 and in section 6, we give conclusion. 

2 Assumptions 

The latest 2003 version of the Wilson program was ap- 
plied for photometric and spectroscopic solutions. The 
method assumes the star surfaces to be equipotential 
and uses a set of curve-dependent or curve- independent 
parameters that can be adjusted by LC and DC pro- 
grams: the orbital inclination ?'. surface potentials O1.2, 
the mean surface effective temperatures Ii^2, the mass 
ratio q — m2/m\, the bandpass luminosities ^1,2, the 
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Table 1 The limb darkening coefficients for EE Aqr 



Parameters 


Filter B 


Filter U 


Filter V 


Xi{bol) 


0.640 


0.640 


0.640 


yi{bol) 


0.255 


0.255 


0.255 


y2{bol) 


0.153 


0.153 


0.153 


Xi 


0.690 


0.795 


0.788 




0.800 


0.802 


0.834 


yi 


0.291 


0.328 


0.279 


2/2 


-0.014 


-0.411 


-0.189 



wavelength-specific limb darkening coefficients x-i_^2, the 
bolomctric limb darkening coefficients xi,2(bol), the 
bolometric gravity darkening exponents gi,2 and the 
bolometric albedos ^1,2- Throughout this paper, the 
subscripts 1 and 2 refer to the primary (hotter) and 
the secondary (cooler) components, respectively. 

For both components of the systems, we used bolo- 
metric linear, logarithmic and square root law and the 
best result was obtained for bolomctric logaritmic limb 
darkening law of Klinglesmith & Sobieski (1970) of the 
form: 

/ = /o(l — X + xcos6 — y cos^ln(cos^)), (1) 

where the limb darkening coefficients x and y for both 

components were fixed to their theoretical values, inter- 
polated using Van Hamme's (1993) formula which have 
tabulated in tables 1 and 2. 

The gravity darkening exponent from Lucy (1967) 
and the bolometric albedos from Rucinski (2001) were 
chosen for convective envelopes (<? = 0.32, A = 0.5) 
and for radiative envelopes {g = 1.0,A= 1.0), which 
are agreement with the final surface temperature. In 
order to reduce the number of free parameters these 
parameters were kept constant during all the iterations. 
Also it is asstimcd that this binary system has zero 
orbital eccentricity (e = 0.0) and that its rotational 
and orbital spins are synchronous (Fi = F2 = 1.0). 
Also black body models arc employed, and we assume 
that there is no third light (^3 = 0.0) for both systems. 

3 Photometric Solutions 

Using the photometric data from the WI and BR, we 

derived the final elements of the systems. Before be- 
ginning the analysis, we have chosen some parameters 
of the systems using the spectroscopic and photometric 
information as a starting input values. 

Initially, the light curve program (LC) was imple- 
mented in mode 2 with no third light (based on the 



Table 2 The limb darkening coefficients for Z Vul 



Parameters 


Filter B 


Filter U 


Filter V 


Xi{bol) 


0.762 


0.762 


0.762 


X2{bol) 


0.703 


0.703 


0.703 


yi{bol) 


0.090 


0.090 


0.090 


y2{bol) 


0.072 


0.072 


0.072 


Xi 


0.434 


0.482 


0.509 


X2 


0.584 


0.615 


0.682 


Vi 


0.233 


0.228 


0.275 


2/2 


0.290 


0.249 


0.336 



spectroscopic observations, HK and PO) correspond- 
ing to detached configuration by choosing i, T2 (while 
the temperature of the primary component was as- 
sumed from spectroscopic or color data), Q.\^2 and Li 
as adjustable parameters in UBV filters. The symme- 
try between the maximums in the light curves indicate 
the lack of star spots on the components. After a few 
runs of the LC program an initial set of values was cm- 
ployed as input parameters for the DC program. After 
a few runs with the differential correction program, we 
could not obtain the values of the physical parameters 
with less error, thus making change to mode 5 so that 
expected from the earlier workers. By several iterations 
of the LC and DC programs and adjusting the param- 
eters i, T2, il\ and Li for each filter, finally the solution 
evolved into a nearcontact configuration for EE Aqr 
and semi-detached one for Z Vul. The results of light 
curve solutions with the final elements are given in Ta- 
bles 3 and 4. The theoretical light curves computed 
with these results arc shown in Figures 1 and 2. The 
agreement between the observed (solid circles) and the 
theoretical light curves (continuous lines) is quite good. 

4 Spectroscopic Solutions 

Using the radial- velocity data from HK and PO, we 
derived the orbital elements of the systems. Adjustable 
parameters were the following: the orbital semi-major 
axis a, the radial velocity of the binary system center 
of mass Ky. The results are given in table 5. Using 
the final elements of the objects, the theoretical radial- 
velocity curves are shown in Figure 3. The agreement 
between thc^ observed (solid circles) and the theoretical 
radial-velocity curves (continuous lines) is quite good. 

5 Absolute Elements 

Using the obtained results of the light and radial- 
velocity curves, we calculate the absolute parameters 



Table 3 The photometric parameters of EE Aqr 

Parameters This work This work This work WiUiamon PadaUa Russo and Hilditch Covino et Covino et 

B Filter U Filter V Filter (1974) (1979) Sollazzo and King al. (1990) al. (1990) 

(Johnson) (Johnson) (Johnson) (1982) (1988) BVRI BVRIkRV 



i 


80.5074 


81.3365 


78.9509 


73.71 


74.00 


76.0 


80.5 


80.0 


80.2 




+ f) 9Q7^ 




















U.OZi 






u.ouo 




n 4 


n ^9 


n ^^9 


n ^97 






















J 1 


1 uuu 






snnn 


snnn 


snnn 


7nfin 


79^n 


7997 




fixed 


fixed 


fixGci 














To 






41 7^ 






4440 




4240 


49^^ 




zb 116 


± 252 
















0,1 




2.6141 


9 (sn^n 

Zi.UUOU 




nn^ 


9 7S 




9 fil Q 






di 0.0118 




di 0.0116 














On 


9 "linn 


9 "linn 




2.882 




9 fiS 




9 ^^S 




r, /(T,. _i_ r^^ 

-1^1/ ^ J-12 ) 






n Q71 

u.y f J-U 


n Qfifi 

U. JUU 




n Q^i 


n Q^9 


n Q^Q 








\J.\J\J0\J 


0984 


n n^4 




n n4Q 


n n4s 


n n4i 








4 


4 














log g2(CGS) 


4.24 


4.24 


4.24 














ri (pole) 


0.4303 


0.4310 


0.4330 


0.394 


0.394 




0.434 


0.43 


0.44 




± 0.0021 


± 0.0019 


± 0.0021 














vi (point^ 


0.5121 


0.5139 


0.5194 


0.469 


0.467 






0.52 


0.54 




± 0.0056 


± 0.0051 


± 0.0059 














ri{side) 


0.4568 


0.4577 


0.4603 


0.415 


0.414 






0.46 


0.47 




± 0.0027 


± 0.0025 


± 0.0027 














nipack) 


0.4765 


0.4775 


0.4806 


0.437 


0.435 






0.48 


0.49 




± 0.0033 


± 0.0030 


± 0.0033 














r2{pole) 


0.2659 


0.2659 


0.2659 


0.300 


0.299 




0.262 


0.27 


0.26 




± 0.0025 


± 0.0017 


± 0.0032 














r2{point) 


0.3853 


0.3853 


0.3853 


0.406 


0.405 






0.36 


0.33 




± 0.0121 


± 0.0080 


± 0.0153 














r2{side) 


0.2769 


0.2769 


0.2769 


0.313 


0.312 






0.28 


0.27 




± 0.0027 


± 0.0018 


± 0.0034 














r2{back) 


0.3096 


0.3096 


0.3096 


0.345 


0.345 






0.31 


0.30 




± 0.0027 


± 0.0018 


± 0.0034 














OiYi 


2.5100 


2.5100 


2.5100 














Oout 


2.3114 


2.3114 


2.3114 














S u;{o - 


0.0172 


0.0234 


0.0156 
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Table 4 The photometric parameters of Z Vul 

Parameters This work This work This work Popper Cester et al. (1977) 

B Filter U Filter V Filter (1957b) B Filter U Filter V Filter 
(Johnson) (Johnson) (Johnson) (Johnson) (Johnson) (Johnson) 



i 


88.6106 


88.5185 


88.6478 88 


88 


88.5 


88.9 




± 0.0666 


± 0.0453 


± 0.0772 








q 


0.43 


0.43 


0.43 


0.43 


0.43 


0.43 




fixed 


fixed 


fixed 








Ti 


19840 


19840 


19840 


19850 


19850 


19840 




fixed 


fixed 


fixed 










10909 


9882 


10810 


9000 


10290 


9410 




± 25 


± 32 


± 27 










3.7975 


3.8512 


3.7434 










± 0.0106 


± 0.0083 


± 0.0119 








^2 


2.7387 


2.7387 


2.7387 








Li/(Li+L2) 


0.8151 


0.8869 


0.7889 


0.962 


0.930 


0.954 


L2/{Li+L2) 


0.1849 


0.1131 


0.2111 


0.038 


0.070 


0.046 


log gi(CGS) 


3.87 


3.89 


3.86 








log g2(CGS) 


3.49 


3.49 


3.49 








ri {pole) 


0.2954 


0.2908 


0.3002 










± 0.0009 


± 0.0007 


± 0.0011 








T\ {point) 


0.3086 


0.3031 


0.3144 










± 0.0011 


± 0.0008 


± 0.0013 








Vi {side) 


0.3011 


0.2962 


0.3063 










± 0.0010 


± 0.0007 


± 0.0012 








ri (back) 


0.3058 


0.3006 


0.3113 










± 0.0011 


± 0.0008 


± 0.0012 








r2{pole) 


0.2881 


0.2881 


0.2881 










± 0.0004 


± 0.0004 


± 0.0005 








r2{point) 


0.4142 


0.4142 


0.4142 










± 0.0018 


± 0.0016 


± 0.0021 








r2{side) 


0.3004 


0.3004 


0.3004 










± 0.0004 


± 0.0004 


± 0.0005 








r2 {back) 


0.3330 


0.3330 


0.3330 










± 0.0004 


± 0.0004 


± 0.0005 








^in 


2.7387 


2.7387 


2.7387 








^out 


2.4781 


2.4781 


2.4781 








S UJ{0 - C)2 


0.0039 


0.0057 


0.0031 
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Table 5 The spectroscopic parameters 



Parameters 


EE Aqr 


Z Vul 


ai sin i = 


V^{Km/s) 


-1.6786 


-22.7690 






± 0.6520 


± 0.3952 




aiRe) 


4.1240 


15.8411 


02 sin i = 




± 0.0188 


± 0.0454 





of the systems. The results are listed in tables 6 and 
7. Wc determine the absolute dimensions using the fol- 
lowing formulae: 



Mi/Mq = /i(Mi, M2, i)(l + (z)7 sin^ i 



/i(Mi, M2, t) = (M2 sinz)V(Mi + M^f 



M2/M0 =g(Mi/M0), 



(2) 



(3) 



(4) 



Ri,2/Re = ^.2m{Mi^2lMQ{l + q)P\days)y'^r{side), 

(5) 



Li,2/i0 = (i^l,2/i^0)'(^l,2/^0)^ 



(6) 



Pi,2/p© = (0.01344Mi,2)/[(Mi + M2)P'rUl (7) 



(Mw)i,2 = 42.36 - 101ogTi,2 - 51og(i?i,2/i?0). (8) 



6 Conclusion 

Comparing the new photometric-spectroscopic solu- 
tions of the systems with those of given in literature, 
the present authors can get following conclusions: 

1-In our study, we may conclude that the systems 
are near-contact and semi-detached for EE Aqr and Z 
Vul respectively. The fillout for the components can be 
calculated from the following formula: 



fillouti^2 



^1: 



X 100. 



(9) 



Also from the following formulae we calculate K12 
and ai,2 sini for each of the systems: 



^ 47r2 J 



ai sm ^ 



(11) 



(12) 



The derived values of the parameters are given in 
tables 8 and 9. According to our results, EE Aqr is a 
near-contact system which the primary and secondary 
components filling are almost 96 ± 1 and 100 percent 
of their respective critical Roche lobes. Also Z Vul 
is a semi-detached system that the primary and sec- 
ondary components filling are almost 72 ± 1 and 100 
percent of their respective critical Roche lobes. The 
fillout percentage was not accurately determined in pre- 
vious workers (86 ± 3 and 94 ± 3 percent according to 
HK for the primary and secondary components of EE 
Aqr, respectively; no information for ones of Z Vul). 
It seems that fillout obtained from our solutions is ap- 
propriate for highly evolved system. 

2- According to our synthetic light curves, EE Aqr is 
(3-Lyrae type and Z Vul is an Algol type system. Also 
according to the obtained parameters of the systems, 
we have drawn the configuration of the components us- 
ing the Binary Maker 2.0 (BradStreet, 1993) software, 
which are shown in figures 4 and 5. 

3- The new values of the velocity amplitudes and 
the gama velocity differ only slightly from earlier re- 
searches. 

4- The absolute elements of both components de- 
termine the evolutionary state of the systems. En- 
tering the results in the mass-luminosity (M-L), Mboi- 
mass (M-M), mass-radius (M-R) and H-R diagrams 
(the diagrams are displayed in figures 6 and 7, it ap- 
pears that EE Aqr contains two main-sequence stars 
with F2IV-I-K4III spectral type while the primary of 
Z Vul is still a main-sequence object, and the sec- 
ondary component is on its way to the giant stage 
with B2V-I-B9V spectral type for the primary and sec- 
ondary respectively (according to tables of Straizys & 
Kuriliene, 1981). The results are tabulated in table 10 
and compared with the results of the pervious authors. 
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Table 6 The absolute elements of the binary system EE Aqr 



ParcLineters 


L Ills vvorK 
J? lllCl Jj 

(Johnson) 


X nis vvorK. 
(Johnson) 


1 niS VVOiK. 
r llXcI V 

(Johnson) 


Russo and 

DOlld-ZzO 
f 1 QSO'^ 


Hilditch cind. 
iviiig ^^lyooj 


Covino et 

dl. ylaiyUJ 


71 /f / ]\/f 


Z. /Di 


Z. iOi 


Z. / Oi 


1 

i.y 





1 /I 
Z.i4 


]\/fr, / A/f ^ 


w.ooo 


w.ooo 


u.ooo 




71 

U. 1 -L 




Li/Lq 


6.55 


6.55 


6.55 


9.16 




7.8 


L2/ Lq 


0.32 


0.25 


0.32 


0.47 




0.33 


Ri/Rq 


1.88 


1.88 


1.88 


1.58 


1.75 


1.79 


R2/RQ 


1.18 


1.18 


1.18 


1.21 


1.07 


1.06 


Pl/PQ 


0.4123 


0.4099 


0.4030 








PiIpq 


0.5922 


0.5922 


0.5922 








Mi{bol) 


2.55 


2.55 


2.53 








M2{bol) 


5.82 


6.08 


5.85 








MM,,M2,i) 


0.0498 


0.0501 


0.0490 









Table 7 The absolute elements of the binary system Z Vul 



Parameters 


This Work 


This Work 


This Work 


Popper (1957b) 


Cester et al. (1977) 




Filter B 


Filter U 


Fihcr V 








(Johnson) 


(Johnson) 


(Johnson) 






M,/Mq 


6.209 


6.209 


6.209 


5.4 


5.4 


M2/M0 


2.670 


2.670 


2.670 


2.3 


2.3 


Li/Lq 


2631.37 


2543.85 


2720.37 


1850 


2818 


L2/LQ 


252.75 


170.21 


243.75 


175 


162 


Ri/Rq 


4.77 


4.69 


4.85 


4.7 


4.5 


R2/RQ 


4.89 


4.89 


4.89 


4.7 


4.6 


Pi/Pq 


0.0571 


0.0600 


0.0543 






P2/P0 


0.0247 


0.0247 


0.0247 






Mi{hol) 


-3.96 


-3.93 


-4.00 


-3.5 




M2{hol) 


-1.42 


-0.99 


-1.38 


-1.0 






0.2412 


0.2412 


0.2412 






Table 8 Other elements of the binary system EE Aqr 






Parameters 


This Work 


This Work 


This Work 


Hilditch and 






Filter B 


Filter U 


Filter V 


King (1988) 






(Johnson) 


(Johnson) 


(Johnson) 






ai sini 


0.9876 


0.9899 


0.9827 


0.92 




02 sin i 


3.0880 


3.0951 


3.0728 


2.85 




Ki{Km/s) 


98.2403 


98.4678 


97.7579 


91.3 




K2{Km/s) 


307.1819 


307.8929 


305.6733 


284 




fillout 1(%) 


95.882 


96.020 


96.428 


86 




fillout 2(%) 


100.000 


100.000 


100.000 


94 
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Table 9 Other elements of the binary system Z Vul 



Parameters 


This Work 
Filter B 
(Johnson) 


This Work 

Filter U 
(Johnson) 


This Work 
Filter V 
(Johnson) 


Plaskett (1920) 


Popper (1957b) 


ai sin i 


4.7716 


4.7714 


4.7717 






02 sin i 


11.0963 


11.0958 


11.0964 






Ki{Km/s) 


98.4167 


98.4127 


98.4182 


96.4 


89.8 


K2{Km/s) 


228.8649 


228.8557 


228.8685 


213.7 


219.7 


fillout 1(%) 


72.118 


71.112 


73.161 






fillout 2(%) 


100.000 


100.000 


100.000 







Table 10 Spectral type of EE Aqr and Z Vul 



Star 


This work 


Williamon 


Padalia 


Russo and 


Plaskett 


Petrie 


Roman 


Poper 


Ccstcr ct al. 






(1974) 


(1979) 


Sollazzo 
(1982) 


(1920) 


(1950) 


(1956) 


(1957b) 


(1977) 


EE Aqr (primary) 


F2 IV 


FO 


F V 


A8 V 












EE Aqr (secondary) 


K4 III 




A 


K3-K4 












Z Vul (primary) 


B2 V 








B3 


B4 


B5 V 


B3-4 V 


B2 V 


Z Vul (secondary) 


B9 V 








B3 


B6 


A 


A2-3 III 


Al III 



Short article title 
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Fig. 1 The observed and theoretical light curves of EE 
Aqr. Solid circles show the observed data and the theoreti- 
cal light curves are shown by continuous lines. 



Fig. 2 The observed and theoretical light curves of Z Vul. 
Solid circles show the observed data and the theoretical light 
curves are shown by continuous lines. 
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Fig. 4 The observed and theoretical radial-velocity curves 
of AV Del. Solid circles show the observed data and the the- 
oretical radial- velocity curve is shown by continuous lines. 
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Phase = 9. 75 

Fig. 6 Configuration of binary system Z Vul 



Fig. 5 Configuration of binary system EE Aqr 
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Fig. 7 M-L and H-R diagram. The location of components 
are shown with star sign. 
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Fig. 8 M-M and M-R diagram. The location of compo- 
nents are shown with star sign. 
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